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^N ■ PACS. 47. 20. Ma - Hydrodynamic stability and instability. 

r^ \ PACS. 75. 50. Mm - Magnetic liquids. 

Abstract. - The dependence of the maximal growth rate of the modes of the Rosensweig 

instability on the properties of the magnetic fluid and the external magnetic induction is studied. 

^AJ ' An expansion and a fit procedure are applied in the appropriate ranges of the supercritical 

Ph , induction B. With increasing B the scaling of the maximal growth rate changes from linear 

r^ ■ to a combination of linear and square-root-like scaling. The scaling of the corresponding wave 

number alternates from quadratic to primarily linear. For very small B the dependence of the 

maximal growth rate on the viscosity is given. Suggestions are made for experiments to test 

the predicted scaling behaviours. 
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(<— ^ . Introduction. - The investigation of instabilities in magnetic fluids has a long history 

[ — ' where the most prominent instability being the normal field or Rosensweig instability ptH]. 
^D . Above a threshold of the induction, the initially flat surface exhibits a stationary array of 
peaks. Despite its long history, some aspects of the Rosensweig instability have been addressed 
only recently: the hexagon-square transition y] or the wave number selection problem ]^^. 

^ ' The wave number is the absolute value of the wave vector, q — |q|, which characterizes small 
disturbances. The ground state of a pattern forming system is subjected to such small distur- 

^ ' bances in order to study its stability. The corresponding growth rate uj may depend on system 

^ . parameters {V}, e.g. the viscosity of the fluid, and control parameters {TZ\, e.g. an external 
magnetic field, lo = Lo{q\{V},{Tl}). Generally it is assumed that in the linear stage of the 
pattern forming process the wave number with the largest growth rate will prevail. Therefore 

C^ ' this mode is called linearly most unstable mode. Due to its role in the pattern formation it 
is of particular interest to examine the maximal growth rate ujrm and its dependence on the 
different parameters. 

This dependence has received only limited attention in classical hydrodynamic systems. 
For the Kiippers-Lortz (KL) instability in Rayleigh-Benard convection rotated about a vertical 
axis, the growth rate of different KL angles were calculated for one fixed rotation rate and 
different temperature differences S. In surface-tension-driven Bcnard convection the growth 
rates are calculated for two fixed values of heat loss and different temperature differences H . 
But in both systems it was not analysed how uom depends on the control parameters. In the 
problem of convection for autocatalytic reaction fronts, the maximal growth rate was analysed 
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only for the two distinct limits of infinite and zero thermal diffusivity M . In both cases a power 
law was found relating the increase of the maximal growth rate with the density difference 
in the fluid. For the Rayleigh- Taylor instability of superposed incompressible viscous fluids, 
the maximal growth rate was analysed for appropriately scaled densities and viscosities in S . 
The dependence of ujm on the remaining parameter, the scaled surface tension, is given by 
a single function covering the whole range of possible values of the surface tension. With 
respect to experiments, no measurements focusing on the maximal growth rate have yet been 
undertaken. 

Particularly in magnetically or electrically driven systems the fastest growing mode is 
experimentally accessible. If the field is increased with a sudden jump, a pattern characterized 
by the linearly most unstable mode should be observed. The validity of this conclusion was 
shown recently for the Rosensweig instability in magnetic fluids Q . The lack of studies and 
the experimental access motivates this letter, where the dependence of the maximal growth 
rate on the properties of the magnetic fluid and the external magnetic induction is studied. 

System. - A horizontally unbounded layer of an incompressible, nonconducting, and 
viscous magnetic fluid of thickness h and constant density p is considered. The fluid is bounded 
from below by the bottom of a container made of a magnetically impermeable material and has 
a free surface with air above. The electrically insulating fluid justifies the stationary form of 
the Maxwell equations, which reduce to the Laplace equation for the magnetic potentials in the 
region of the container, the magnetic fluid, and the air. It is assumed that the magnetization 
of the magnetic fluid depends linearly on the applied magnetic field, M = (/.t^ — 1)H, where 
Hr is the relative permeability of the fluid. 

In a linear stability analysis, all small disturbances from the basic state are decomposed 
into normal modes, i.e., into components of the form exp[~i{ujt — qr)] with r = {x,y). If 
Im(w) > 0, initially small undulations will grow exponentially and the originally horizontal 
surface is unstable. Therefore uj is commonly called the growth rate, which is in fact true only 
for its imaginary part in the chosen normal mode ansatz. The linear stability analysis leads 
to the dispersion relation [pl-pT[ 



^ - ,W-M ;wTT 9 [4g^ + iq' + fr] cotHqh) - q [Aq^^ + {q^ + f)'] tanh(g/i) 

(7 coth((7/i) — (7 coth(g/i) ' 



■ , ., . ■ ,,,. > + tanh{qh) 

cosh[qn) sinh[qn) J 



gq 



q'n 


cot'h{qh) — q Aq^if 


-q^ 
P 


fJ-oPvM^ 2" 

Hqh)q 

P 



(1) 



where v is the kinematic viscosity, a the surface tension of the magnetic fluid with air, M 
the absolute value of the magnetization, g = (0, 0, — g) the acceleration due to gravity, po the 
permeability of free space, g = y^q^'^^HoJ^, and A{qh) = [e''''(l + /Zr) +e~''^(l-/Xr)]/[e''''(H- 
IJ,r)^ — e~^''(l — P-r)'^]- For not too thin layers the approximation with an infinitely thick layer 
is applicable MM- Therefore the starting point is the dispersion relation (nh for h ^ 00 |13] 
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Dimensionless quantities are introduced for all lengths, the induction, the time, and the vis- 
cosity, 

I = qd, ^ = 7r-' (3) 
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where tc is the so-caUed capillary time. For w = iW2 with uj2 > 0, the real part of eq. 
reduces to 



2 ;7_i_^3 _ or2;=2 



/+(,-,^,;z.,S):=(. + ^ +^—^-r^f^^~^\l + i^=0. (5) 






The imaginary part of eq. (m) is identically zero. The wave number of maximal growth, g™, 
is defined by doj2/dq — 0. Since u)2 is given implicitly by /+, the cross section of /+ = and 
df-\-/dq — determines the maximal growth rate and its corresponding wave number. 

Analysis and Results. - An expansion of _B, q, and uj2 in the following form 

B =1 + B q = 1 + qrn ^2=0 + Lb2,rn (6) 

leads to an analytical expression of the dependence of t02,m and qm on the induction and the 
viscosity. All hatted quantities in (g) are small, {B,qrmi02,m) ^ 1, and denote the scaled 
distances from the critical values at the onset of the instability. If i? 3> 0J2,m, the expansion 
of /+ = and its derivative results in 

A9uJ2,.n - 8B + 3loI^+ 8iyq„,L02,rn - l6qn.B - AB^ - ^ = (7) 

^3 

- 16B - 8^2 + 4(7™ - 16Bg™ + 8^(2.2,™ + 8P<7„(i2,m + ^^ =0. (8) 

Considering linear terms in eqs. (R) and (p|) only, one obtains 0)2, m = (2/i^)i3 and qm = 0. The 
latter equation describes an incorrect dependence of the wave number of maximal growth on 
the applied induction. The result g,„ = and q = 1 (sec eq. (I^)), respectively, implies that 
the wave number of maximal growth is constant and equal 1 with increasing induction. Such 
a conclusion is true only in the case of infinitely viscous magnetic fluids as shown in |j] . Now 
higher order terms of the applied induction arc included by the ansatz 

LU2^„, ^ aB + I3B^ + 7^3 + 0{B^) , (9) 

q,n = SB^ + (B^ + 0(5^) . (10) 

Using such an ansatz, the equations (R) and (@) contain all terms that contribute up to 
third order in B. By determining the expansion coefficients a, . . . , e, the dependence on the 
parameters viscosity and induction is now given by 

2-/1 3\- /10 3\- 
W2,m = -B +[-- — ] B"^ +[ — - — ] B"^ + 0{B^) for Q < B < d'^/Q (11) 

qm - 4^' +(^~^)b' + 0{B^) for Q<B< pV6 ■ (12) 

For scaled inductions larger than P^ /6, one has to solve the full implicit equation (|l|) and its 
derivative with respect to q numerically. Through the implicit character of both equations, 
parameter fits are possible only for the dependence of 6j2,m and qm on B. The fit for an 
excellent agreement with the numerical data includes a linear term and a square-root term 
with respect to B, where the coefficients depend on the used magnetic fluid. Combining both 
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Table I - Material ||13[ and fit parameter of the magnetic fluids. 



Fluid 


P (kg 


EMG 909 


1.02 


EMG 901 


1.53 


APG J16 


1.01 


APG J14 


1.06 


APG S20 


1.05 


APG J12 


1.11 


APG L17 


1.05 


APG JIO 


1.16 


APG Slln 


1.15 


APG 077n 


1.19 


APG 037 


0.92 



lo-" 

10^ 
10^ 
10^ 
10^ 
10^ 
10^ 
10^ 
10^ 
10^ 
10^ 



5.880 
6.540 
2.475 
2.830 
3.333 
3.604 
5.714 
6.034 
6.957 
9.244 
1.196 



-') ^ (kg s-2) 



10'" 
10"'' 
10-^ 
10"^ 
10"^ 
10-^ 
10"^ 
10-^ 
lO"'' 
10"^ 
10""* 



2.65 
2.27 
3.30 
3.40 
3.30 
3.40 
3.40 
3.40 
3.30 
3.50 
3.40 



10"^ 
10"^ 
10-^ 
10"^ 
10"^ 
10"^ 
10-2 
10-2 
10-2 

10^2 
10^2 



/?2/6 
1.4- 10-^ 

• 10-" 

■ 10"^ 

■ 10"^ 
•10"^ 

• 10"^ 

■ 10-^ 

■ 10-^ 
■10-3 

■ 10^3 

■ 10^2 



3.9 
1.7 
2.3 
3.3 
4.0 
9.2 
1.2 
1.6 
2.8 
3.3 



Cl 

1.03 
1.16 
0.72 
0.69 
0.74 
0.58 
0.43 
0.38 
0.33 
0.26 
0.25 



C2 

3.90 
2.95 
3.52 
3.44 
2.91 
3.29 
3.02 
2.92 
2.80 
2.55 
2.43 



3.46 
3.32 
3.12 
3.08 
2.99 
2.95 
2.71 
2.63 
2.51 
2.34 
2.24 



C4 

-0.07 
-0.10 
-0.20 
-0.20 
-0.22 
-0.23 
-0.26 
-0.27 
-0.28 
-0.32 
-0.32 



ranges oi B, we finally have 
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for < B < pV6 

for j?V6 < S < 0.4 

for < B < ^2/6 

for i>V6< S <0.4, 



(13) 
(14) 
(15) 
(16) 



where the four coefficients Ci, i = 1, . . . ,4, are given in table | for eleven different fluids. These 
magnetic fluids are made of magnetite nanoparticles dispersed in a carrier liquid (synthetic 
ester for APG fluids, light petrol for EMG fluids). The nanoparticles are coated with a layer 
of chemically adsorbed surfactants to avoid agglomeration. 

From eqs. (|l^l6) and the figures ^^ it is evidently that there are two different scaling 
regimes for the maximal growth rate as well as for the corresponding wave number. In all three 
examples, the assumption for the expansion, v ^ 0)2, m, is well accomplished, see horizontal 
arrows and solid lines in figs. |l|(a)-|^(a). A lower bound for the transition between the two 
regimes is given by P^/G. This value is denoted by the vertical arrows in figs. yj3|. i/ is a 
combination of the density, the viscosity, and the surface tension of the magnetic fluid, see 
eq. (y). Therefore the transition range can be matched to experimentally manageable step sizes 
for the jumplikc increase of the magnetic induction by choosing a proper fluid. For the chosen 
three fluids, the transition range shrinks from a broad interval of 0.0005 < Btrans < 0.004 
for EMG 901 (fig. §) to Btrans ~ 0.03 for APG 037 (fig. |). This shift is mainly caused 
by an increase in the viscosity (see table ffl). As a consequence, the range of validity of the 
two-parameter fit in eqs. (14) and (16) can be replaced by D^/6 < B < 0.4 for the fluid APG 
037 (see fig. H). The appearance of different scaling regimes in the Rosensweig instability 
make it more attractive than instabilities with a single scaling such as the Rayleigh- Taylor 
instability |^. 

The maximal growth rate starts to increase linearly with B. Towards the upper bound 
of the expansion region, B < i^^ /6, corrections appear due to the quadratic and cubic term 
in (ttSh. For D^/6 <^ B, the dependence of the maximal growth rate on the induction is 
given by a linear as well as by a square-root term. Both terms are equally important since 
approximations with only one term result in inferior fits. 
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For very small B, the wave number with maximal growth rate increases quadratically with 
B. Near the upper bound of the expansion region cubic corrections occur. For i/^/6 ^ B, q^ 
depends primarily linearly on B because C3 3> |c4|. This linear dependence was already studied 
in W and a good agreement was found between the measured values and the theoretical data. 

The fits (14) and (16) give only the dependence of uj2.m and qm on the scaled magnetic 
induction. It would beneficial to gain insight into the dependence on i/, too, at least approx- 
imately. Therefore the set of eleven fluids and their fit coefficients is used to approximate 
the relations Ci = Ci(P). It is apparent from fig. H that all coefficients exhibit a smooth de- 
pendence on D for 0.03 ^ i^ S 0.45. Only for D ~ 0.04 and ~ 0.14, the values of ci and 
C2 show deviations. Since 0)2. ,„ ^ <^^ ^^2,m-, inaccuracies in the surface tension a may result 
in incorrect values of ci and 02- From experiments with surface wave damping on ordinary 
fluids it is known that deviations in the value of the surface tension due to contamination 
may account for differences in the order of 20% [^;|l^. Excluding the large deviations, the 
following dependencies of Ci on D are suggested (see solid lines in fig. ^ 



Cl = 


= 1.11- 


- 4:.13D + 5.09^ 


C2 = 


= 4.0 - 5.09 + 3.449^ , 


(17) 


C3 = 


= 3.53- 


-4.219 + 3.069'^ 


C4 = 


= -0.75 + 0.38/^. 


(18) 



It is emphasized that the data cover a range of about 2 orders of magnitude in the kinematic 
viscosity i^. 

Comparing the numerical data and the analytical results for /i ^ 00 of g^, one realizes 
that with increasing fluid viscosity thicker layers are necessary for a good agreement. This 
is illustrated at a fixed value oi B = 0.0002 for the three different fluids. For EMG 901 the 
analytical value according to eq. (15) is already reached at a flnite thickness of /i = 7 mm. 
For APG S20 the layer thickness has to be 14 mm before matching with the value of (jm for 
h -^ 00. Finally, for APG 037 a thickness of ft, = 19 mm is needed to reach qm in the inflnitely 
thick case. Such a sensitivity does not occur for the maximal growth rate. 

Conclusion. - The dependence of the maximal growth rate on the properties of the mag- 
netic fluid and the external magnetic induction is studied. An expansion and a fit procedure 
are applied in the appropriate ranges of the supercritical induction. Both the maximal growth 
rate LU2.m and the corresponding wave number q^ follow different scaling regimes. Depending 
on the material parameters the transition range and its width can be chosen. The scaling of 
'jJ2,m with respect to B changes from linear to a combination of linear and square-root-like 
with increasing B. The scaling of q^ with respect to B alternates from quadratic to essentially 
linear with increasing B. For very small B, Cj2.m depends inversely on the viscosity v of the 
fluid, whereas qm varies as v^^ with the viscosity, see eqs. (^, |l^, 15). 

The use of magnetic fluids with properties as APG S20 are suggested to test experimen- 
tally the predicted scaling regimes. The range, where the expansion holds lies within the 
experimentally realizable step size for the jumplike increase of the induction. A clear sep- 
aration between the different scaling regimes is present (see fig. 0). Thus each regime can 
be examined independently. Since magnetic fluids are opaque and feature poor reflectivity, a 
radioscopic detection technique is the best method for this purpose. By utilizing the attenua- 
tion of X-rays, the dynamic and static properties of the surface deformation can be measured 
with high accuracy jlq]. Applying this method to the problem presented here would fill the 
gap of experiments measuring especially maximal growth rates. 

* * * 
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Fig. 1 - Scaled maximal growth rate 0)2, m (a) and scaled maximal wave number q™ (b) as a function 
of the scaled supercritical induction B for EMG 901. The solid lines denote the analytical results 
(O) for Cj2.m and (13) for ijm in (a) and (b). The long-dashed lines denote the fits (12) for 0)2, m and 
(14) for (jm in (a) and (b). Both scaling regimes are clearly separated. The data are calculated for 
h = 100 mm (o), 50 mm (*), 10 mm (+), and 4 mm (n) from eq. (|l|). The vertical arrows in (a) and 
(b) indicate t'^/G, the horizontal arrow in (a) the dimensionless viscosity v. 
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Fig. 2 - Scaled maximal growth rate a>2.m (a) and scaled maximal wave number Qm (b) as a function 
of the scaled supercritical induction B for APG S20. The transition range between the different 
scaling regimes shifts to higher values of 13 and its width shrinks. The symbols, the arrows and the 
line types are used as in fig. hi. 
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Fig. 3 - Scaled maximal growth rate il>2,m (a) and scaled maximal wave number g™ (b) as a function 
of the scaled supercritical induction B for APG 037. The transition range between the different 
scaling regimes shifts to even higher values of B. Both scaling regimes overlap around B ~ 0.03. The 
symbols, the arrows and the line types are used as in fig. hi. 
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Fig. 4 - Dependence of the coefficients ci (a), C2 (b), cs, (c), and C4 (d) on tlie scaled viscosity v. Tfie 
values of Ci from table nI are plotted as filled squares. The solid lines indicate the approximations by 
eqs. (|l|) and (|l|). 



